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ABSTRACT 

We present Spitzer observations of the dark globule CB190 (L771). We observe a roughly circular 
24/im shadow with a 70" radius. The extinction profile of this shadow matches the profile derived 
from 2MASS photometry at the outer edges of the globule and reaches a maximum of ^ 32 visual 
magnitudes at the center. The corresponding mass of CB190 is ^ 10 Mq. Our ^^CO and ^"^00 J 
= 2-1 data over a lO'xlO' region centered on the shadow show a temperature ^ 10 K. The thermal 
continuum indicates a similar temperature for the dust. The molecular data also show evidence of 
freezeout onto dust grains. We estimate a distance to CB190 of 400 pc using the spectroscopic parallax 
of a star associated with the globule. Bonnor-Ebert fits to the density profile, in conjunction with this 
distance, yield £,max — 7.2, indicating that CB190 may be unstable. The high temperature (56 K) of 
the best fit Bonnor-Ebert model is in contradiction with the CO and thermal continuum data, leading 
to the conclusion that the thermal pressure is not enough to prevent free-fall collapse. We also find 
that the turbulence in the cloud is inadequate to support it. However, the cloud may be supported 
by the magnetic field, if this field is at the average level for dark globules. Since the magnetic field 
will eventually leak out through ambipolar diffusion, it is likely that CB190 is collapsing or in a late 
pre-collapse stage. 

Subject headings: ISM: globules - ISM: individual (CB190) - infrared: ISM - (ISM:) dust, extinction 



1. INTRODUCTION 

Cold cloud cores, where star formation begins, repre- 
sent the stage in early stellar evolution after the forma- 
tion of molecular clouds and before the formation of Class 
objects. Their emission is inaccessible at shorter wave- 
lengths, such as the near infrared and visual bands, due 
to low temperatures, very high gas densities, and associ- 
ated large amounts of dust. Because cold cloud cores can 
best be observed at sub- millimeter and far-infrared wave- 
lengths, these spectral regions are essential to developing 
an understanding of the fir st step s toward star-formation 
(see, e.g.. iBacmann et all l2000t IxTrk et al.ll2007f) . The 
wavelength range accessible to the Spitzer Space Tele- 
scope, 3.6 /im to 160 /im, is ideally suited to observe 
cold, dense regions. 

CB190 (L771) is an example of one such dark globule 
and i s classified in the Lynds catalo g as having an opacity 
of 6 (lLvndslll962l ^. i.e., very high. IClemens fc BarvainisI 
(|1988( ) studv this object as part of an optically se- 
lected survey of small molecular clouds. They find that 
it appears optically isolated, is somewhat asymmetric 
(a/6 ~ 2.5) and has some bright rims of reflection and 
Ha. CB19 is ~ 5^ across, an d has an estimated distance 
of 400 pc (jNeckel et al.l[l980l) . 

We present Spitzer maps of CB190. In particular, we 
highlight the observation of this globule in absorption 
at 24 //m. We combine these data with SCUBA obser- 

^ This work is based in part on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, CaUfornia Institute of Technology, under NASA con- 
tract 1407. 

^ Department of Astronomy and Steward Observatory, Univer- 
sity of Arizona, 933 North Cherry Avenue, Tucson, Arizona 85721; 
astutz@as.arizona.edu. 

^ Jet Propulsion Lab, California Institute of Technology, 4800 
Oak Grove Drive, Pasadena, CA 91109. 



vations at 850 /im. We have obtained complementary 
Hcinrich Hertz Telescope (HHT) ^^CO and "CO J=2-l 
on the fly (OTF) maps of this globule and have used the 
HHT and Green Bank Telescope (GBT) to measure high 
resolution line profiles for ^^CO, ^^CO, NH3, CCS, C3S, 
and HC5N. We report C^^O and DCO"*" measurements 
with the Caltech Submillimeter Observatory (CSO). We 
discuss the issue of stability and possible support mech- 
anisms in some detail; however, we cannot say conclu- 
sively if CB190 is in equilibrium. In § 2 we describe the 
observations and data processing. In § 3 we present our 
main analysis of CB190: we derive an optical depth and 
an extinction profile for the 24 /im shadow using a tech- 
nique presented for the first time in this work'*. We also 
discuss two sources associated with CB190 and derive a 
distance estimate. In § 4 we compare various mass esti- 
mates of this object. In § 5 we describe our Bonnor-Ebert 
fitting method and discuss possible support mechanisms 
for CB190. Finally, in § 6 we summarize our main con- 
clusions. All positions are given in the J2000 system. 

2. OBSERVATIONS AND PROCESSING 
2.1. Spitzer data 

Object CB190 (L771), centered at about RA = 
19'*20"48", Dec = -|-2 3°29^45" was observed with the 
MIPS instrument f.Rieke et al.l 12004) at 24 /im, 70 /tm 
and 160 /fm, Spitzer program ID 53 (P.I. G. Rieke). The 
observations were carried out in scan map mode. Fig- 
ure[T]shows these data, along with the Digital Sky Survey 
(DSS) red plate image of the globule. 

The 24 /im data were reduced using version 3.06 of the 
MIPS Data Analysis Tool (DAT; Gordon et al. 2005). 
In addition to the standard processing, several other ad- 

* Previo us related work h as been conducted using the ISO 7/im 
band, (e.g.. IBacmann et al.|[20 00) 
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ditional steps have been applied: 1) correcting for vari- 
able offsets among the four readouts; 2) applying a scan- 
mirror position-dependent flat field; 3) applying a scan- 
mirror position-independent flat field to remove long 
term gain changes due to previous saturating sources; 
and 4) background subtraction. For the last two steps, 
masks of any bright sources and also of the region of in- 
terest were used to ensure that the two corrections were 
unbiased. The background subtraction was performed 
by fitting a low order polynomial to each scan leg of the 
masked data and subtracting the resulting fit. This pro- 
cedure removes the contribution of zodiacal and other 
background light as well as small transients seen after a 
boost frame. 

The 70 fim and 160 ^m data were also reduced using 
the DAT (Gordon et al. 2005). After completing the 
standard reduction processing, an additional correction 
was performed to remove the long term drift in the Ge:Ga 
detectors. The correction was determined by fitting a 
low order polynomial to the masked version of the entire 
dataset for each pixel. By masking bright sources as 
well as the region of interest we ensure that the masked- 
version fits are unbiased. The resulting fits for each pixel 
were subtracted to remove the long term drift as well as 
any background light. 

CB190 was ob served with the IRAC instrument 
(jFazio et al.ll2004l ) at 3.6, 4.5, 5.8 and 8.0 ^m, program 
ID 94 (P.I. C. Lawrence), see fig. [5] Standard pack- 
ages were used to reduce the data, and the mosaicked 
frames were generated with the MOPEX software pack- 
age. The data were taken in high dynamic range mode; 
the 30 second exposure time was divided into a 1.0s 
"short" frame and a 26.8 s "long" frame at each posi- 
tion. Each observation was repeated 5 times, yielding 
an effective long-frame exposure time of 134 s. SExtrac- 
tor (Bcrtin & Arnouts 1996) was used for both source 
extraction and photometry. The photometry was cross- 
checked with PhotVis version 1.1, an IDL GUI-base d 
implementation of DAOPHOT (|Gutermuth et al.ll2004l ). 
We found good agreement between the two sets of pho- 
tometry. 

2.2. 12(70 and ^^CO data 

The CB190 region was mapped in the J=2-l transi- 
tions of i^CO and "CO with the lO-m diameter HHT 
on Mt. Graham, Arizona on 2005 June 9. The receiver 
was a dual polarization SIS mixer system operating in 
double-sideband mode with a 4 - 6 GHz IF band. The 
12C0 J=2-l hue at 230.538 GHz was placed in the up- 
per sideband and the "CO J=2-l line at 220.399 GHz 
in the lower sideband, with a small offset in frequency 
to ensure that the two lines were adequately separated 
in the IF band. The spectrometers, one for each of the 
two polarizations, were filter banks with 1024 channels 
of 1 MHz width and separation. At the observing fre- 
quencies, the spectral resolution was 1.3 km and the 
angular resolution of the telescope was 32" (FWHM). 

A 10' X 10' field centered at RA = 19''20™49.5^ Dec = 
+23°29'57" was mapped with on-the-fly (OTP) scanning 
in RA at 10" sec"^, with row spacing of 10" in declina- 
tion, over a total of 60 rows. This field was observed 
twice, each time requiring about 100 minutes of elapsed 
time. System temperatures were calibra ted by the stan- 
dard ambient temperature load method (|Kutner fc UlichI 



[1981 after every other row of the map grid. Atmospheric 
conditions were clear and stable, and the system temper- 
atures were nearly constant at Tgys = 450 K (SSB). 

Data for each polarization and CO isotopomer were 
processed with the CLASS reduction package (from the 
University of Grenoble Astrophysics Group), by remov- 
ing a linear baseline and convolving the data to a square 
grid with 16" grid spacing (equal to one-half the telescope 
beamwidth). The intensity scales for the two polariza- 
tions were determined from observations of DR21(0H) 
made just before the OTP maps. The gridded spectral 
data cubes were proce ssed with the Miriad software pack- 
age (jSault et al.lll995l ) for further analysis. The two po- 
larizations were averaged, yielding images with rms noise 
per pixel and per velocity channel of 0.15 K-T^ for both 
the "CO and "CO transitions. 

The linewidths were narrow, so that 70% of the flux in 
the 12 CO J=2-l line was in a single 1 MHz spectrometer 
channel, while essentially all the flux of the ^'^CO line was 
in a single channel. We can therefore set an upper limit of 
~ 1.3 km s^i on the linewidth for the emission lines, but 
have little or no kinematic information from the maps, 
other than the LSR velocity, which i s 11.0 km s~i, in 
agreement with IClemens &: BarvainisI ()1988f ) . In figure [3] 
we show two maps of the integrated ^^CO and ^^CO J=2- 
1 lines, summed over the 2 spectrometer channels with 
detectable emission. Furthermore, in figure [4l we show 
the I'^CO J=2-l contours overlayed on the 24 y^m image; 
the spatial coincidence between the two is evident. 

To better constrain the CO line properties, we observed 
the core of the molecular cloud with high velocity reso- 
lution on 2006 June 21 with the HHT. The position ob- 
served was at RA = 19''20™46.4", Dec = -h23°29'45"6, 
which is the peak of the ^^CO intensity map (fig. [3]). Both 
the 12 CO and ^^CO J=2-l transitions were observed for 
15 minutes each. The high resolution ^^CO and ^^CO 
spectra are shown in fig. [51 The ^^CO line is slightly 
asymmetric, with a peak beam-averaged brightness tem- 
perature of 7.65 K, and a velocity width (FWHM) of 
1.22 km s-i. The "CO hue is narrower, with a FWHM 
of 0.97 km s^i, and has a peak intensity of 2.96 K. The 
i^CO/i^CO intensity ratio at the line peak is therefore 
~2.6; if the isotopic ratio ["CO/"CO] = 50, the optical 
depth of the ^^CO line at the peak is ~0.5, and the line 
is optically thin. While values for ["C0/ "C01 in the 
range of 50 to 70 are reasonable, see e.g., iMilam et al.l 
(|2005f ). changing the isotopic ratio will not significantly 
affect our calculated optical depth. We can then com- 
pute an integrated CO column density assuming the CO 
rotational levels are in PTE. Following [Rohlfs fc Wilswll 
(2004), the peak ^^GO line brightness temperature im- 
plies a CO excitation temperature of 12.6 K. Assuming 
this applies to both isotopomers, the integrated ^^CO 
J=2-l line intensity gives an integrated column density of 
N("CO)= 2x 10" cm-2 at the cloud peak, and N("CO) 
= 9 X IQi^ cm~2. If the [CO/H2] abundance ratio were 
1.5 X 10"'', typical of molecular clouds, the column den- 
sity of H2 would be N(H2) = 6 x 10^° cm'^. A standard 
gas to dust ratio and extinction law would then imply 
Av — 0.7 mag through the cloud core. This very low im- 
plied extinction is clearly incompatible with the observed 
large extinction evident in the POSS image from which 
the L771 dark cloud was identified. This discrepancy 
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suggests that the CO molecule is substantially depleted 
by freeze-out onto dust grains in the core of the cloud, as 
has been seen in many other molecular cloud cores (e.g., 
Tafalla et al. 2002). We confirm this conclusion in § 4.3. 

2.3. Observations of other molecular lines 

Observations of CB190 were performed with the 105 
meter Green Bank Telescope on September 20, 2006. 
Five spectral lines were observed simultaneously in dual 
polarization: NH3 (1,1) and (2,2), CCS Nj = I2 - 2i, 
C3S J = 4 - 3, and HC5N J = 9 - 8. The correlator was 
set up with 6.1 kHz resolution and 8 spectral windows 
(4i?i? and 4LL polarizations) with a 50 MHz bandpass. 
The peak ^^CO position was observed for 20 minutes 
of ON-source integration time while frequency switching 
with a frequency throw of 4.15 MHz. The atmospheric 
optical depth at 1.3 cm was calibrated using the weather 
model of Ron Maddalena (private communication, 2006). 
The average opacity was ri.3 = 0.095 ± 0.010 during the 
CB190 observations. The main beam efficiency was de- 
termined from observations of the quasars 3C286 and 
3C48 and was rjmb = 0.75 ± 0.04. The observations were 
reduced using standard GBTIDL script for frequency- 
switched observations, calibrated using the latest K-band 
receiver Tcai for each polarization, and corrected for at- 
mospheric opacity and the main beam efficiency. 

Michael M. Dunham (private communication, 2006) 
provided C^^O J = 2 - 1 and DCO+ J = 3 - 2 observations 
obtained with the Caltech Submillimeter Observatory. 
The observations were performed in position-switching 
mode with the 50 MHz AOS backend. The main beam 
efficiency was measured to be 0.74 during the observing 
run. C^*0 J = 2 - 1 was detected toward the 24 //m 
shadow peak position, but DCO"*" was not detected to a 
3 a rms of 0.4 K. 

The molecular line emission detected toward CB190 is 
striking in its lack of diversity. Only CO isotopomers 
and NH3 have been detected to date. The early-time 
molecules CCS, C3S, and HC5N were not detected with 
the GBT to a 25 mK (T^) baseline rms while the deuter- 
ation tracer, DCO+ was not detected in the CSO ob- 
servations. The NH3 (1,1) line is weak, with a peak 
T^ = 200 ± 23 mK and a narrow linewidth of Av — 
0.47±0.08 km s^^. Since the main line is a blend of many 
hyperfine lines, the actual linewidth is smaller. The opti- 
cal depth is low enough that the satellite lines are barely 
detected at the 3<t level. Assuming an excitation tem- 
perature of 10 K and optically thin emission, the column 
density of NH3 is TV = 5.4l:*;| x lO^^ ^^-2^ rp^^.^ 
umn density is two orders of magnitude be low the me- 
dian column density in the NH3 survey of iJijina et al.l 
(|l999f l and is six times lower than their weakest detec- 
tion. A modest 20 minute integration time on the GBT 
can probe very low column densities of NH3. Unfortu- 
nately, the (2,2) line was also not detected to the 23 mK 
baseline rms level; therefore, we cannot independently 
determine the kinetic temperature of the gas. The ex- 
tremely weak NH3 emission and DCO+ non-detection 
may indicate that CB190 is a relatively young core. In 
contrast, there is evidence that CO is depleted indicat- 
ing that the core is not a nascent dense core. A more 
extensive and sensitive search for molecular line emis- 
sion should be attempted toward CB190 to characterize 
its chemistry. 



2.4. 2MASS Data 

We used the 2MASS All- Sky Point Source Catalog 
(PSC) (jSkrutskie et al.ll2006f l photometry. We quote the 
default J, H, and K-band photometry in this work, la- 
beled j-m, h_m, and k_m in the 2MASS table header. 
The magnitudes are derived over a 4" radius aper- 
ture. We use the combined, or total, photometric uncer- 
tainties for the default magnitudes, labeled j_msigcom, 
h„msigcom, and k_msigcom in the 2MASS table header. 

2.5. SCUBA Data 

We include in this work the lVisser et al.l ()2002D reduced 
850/xm SCUBA map of CB190 (Claire Chandler, private 
communication, 2006). This map was convolved with 
a 32" FWHM Gaussian beam and is shown in figure [6] 
as contours overlayed on the 160 /im data. The spatial 
agreement between the two wavelengths is very good. We 
note, however, that the southern edge of the cloud may 
be artificially sharpened at 850/im due to the position 
angle of chopping during the scan map. Figure [7| shows 
photometry for the 70 /^m, 160 /im and 850 //m observa- 
tions, measured with a 48" radius aperture centered on 
the 24 /im shadow coordinates. The observed ratio of the 
long wavelength fluxes is /,y[160/.tm]//i,[850^m] = 3.1. 
With this ratio we fit for a cloud temperature using the 
model 

cx vl^B,{T). (1) 

We fix the value of /3 at 1.5 and 2.0, a reas onable range 
of values for dust emissivity ()Whittetll2003f ). and derive 
best-fit model temperatures of 12.0 K and 10.4 K, re- 
spectively. These two models are plotted in figure [71 
The 160 /im flux is well detected, with a signal to noise 
ratio ~20, indicating that these temperatures are robust. 
However, to test the models conservatively, we recalcu- 
late model temperatures allowing for 20% errors in our 
photometry. We do not find significant variation in the 
derived temperatures. We note that the lack of a detec- 
tion at 70 /im mildly favors the colder temperature of 
10.4 K. The 3cr upper limit to the flux density at 70 fxm 
is 16 mjy while the predicted values are 17 mJy for 
the (3 = 1.5, T = 12 K model, and 6 mJy for the (3 = 2.0, 
T = 10.4 K m o del 

iVisser et al.l ([200I flnd that the CB190 SCUBA 
450 /im emission is spatially offset from the 850 /im data 
for reasons that are not understood. Because we flnd 
good spatial agreement between the 850 /im data and 
the 160 /im image (see flg. jS]), we do not use the 450 /im 
map. The need for deep sub- mm or mm observations 
of this region is highlighted by the fact that both of the 
SCUBA maps have low S/N, suffer from ambiguities in 
the spatial extent of the cloud due to the chopping posi- 
tion, and have poorly understood spatial disagreements 
between the 850 /im and 450 /im observations. 

2.6. Optical Spectrum 

We have obtained optical spectra of three stars: 
HD344204 (star 1 in figure H]), HD1608, and Vega. These 
spectra cover the optical range, from 3615A to 6900A 
with an effective resolution of AA = 9A, and were ob- 
served in July 2006 at the Steward Observatory Bok 
2.3m telescope at Kitt Peak using the Boiler and Chivens 
Spectrograph with a 400/mm grating in first order. They 
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were processed with standard IRAF data reduction pack- 
ages. These data are shown in figure [8] and discussed in 
SS 3.1 and 3.2. 

3. ANALYSIS 
3.1. Two associated stellar sources 

There are two bright sources in the 24/j,m image that 
are hkely associated with CB190. The first, labeled 
source 1 in fig.lH is a bright point source, just to the north 
of the 24 /im shadow, with a large amount of diffuse emis- 
sion. This star is HD344204 (IRAS 19186-1-2325), and is 
located at RA = 19''20™47^ Dec = 23°31'40.6". The 
second star, labeled source 2 in fig. IH is spatially coin- 
cident with a small peak in the ^^CO emission, and is 
located at RA = 19''20™57^ Dec = 23°31'37.6". The 
broad-band SEDs of these two sou rces are plot ted in 
fig. El and include 2MASS (Skruts kie et al.l[200l J, H, 
and K data (see § 2.3), IRAC [3.6 /xm], [4.5 /xm], [5.8 ^im], 
and [8.0 /zm] data, and MIPS 24 fj.m fluxes. Although 
star 1 appears to be indistinguishable from the surround- 
ing diffuse emission in figure 21 this is only due to the 
scale used to display the image. To measure the flux 
of the star while avoiding contamination from the sur- 
rounding diffuse emission we use a very small photom- 
etry aperture radius, 6.23", and sky annulus, 6.23" to 
7.47". We derive an aperture correction of 2.1 to this 
flux using an isolated point source, measured with the 
same aperture geometry as that listed above. We com- 
pare this result to the flux derived using a 13" aperture, 
a 20" to 32" sky annulus, and the 24 fim aperture cor- 
rection recommended by the Spitzer Science Center. For 
comparison fig. [§] also indicates the 24 /im fiux of Star 1 
measured with a large aperture with radius = 30" and 
an inner and outer sky annulus radius = 34" and 38" re- 
spectively, meant to include all the light from the diffuse 
emission surrounding the source. 

In fig. |9l we also show the broad-band SED of source 
2. The IRAC colors are ([3.6 //m]-[4.5 ^im\) = 0.5 mag, 
and ([5.8 Atni]- [ 8.0 ii m]) = 1.1 mag. The models by 
iWhitnev et al.l (|2003f ) suggest that these colors are con- 
sistent with those of a late Class source. 

3.2. The Distance 

To understand the nature of the CB190 24 /im shadow 
one must measure its physical properties, such as mass 
and size; to do so one needs an accurate distance. While 
star counts are commonly used to estimate distances to 
nearby clouds, in the case of CB190 this method is not 
reliable due to the small number of foreground sources. 
Another distance estimator is the LSR velocity relation. 
However, this method is not reliable for nearby objects 
whose mot ions are still locally dominated, as is the case 
of CB190. iNeckel et all (|1980[ ) estimate the distance to 
this cloud to be ^400 pc using th e discontinuity in Ay 
with distance. iDame fc Thaddeui ((1985) argue that this 
distance is consistent with the narrow line width they 
measure and therefore reject the other plausible distance 
to this cloud, that of the Vul OBI association at 2.3 kpc, 
noting that this longer distance would imply a much 
larger mass and line width. 

We can obtain a rough estimate of the distance 
using the col o rs of source 2 (see fig. [S]) and the 
IWhitnev etld] (|2003f) models. Their color -magnitude 



([5.8 /im]-[8.0 /im]) vs. [3.6 /im] relation for a face-on late 
Class source, with [3.6 /im] ~ 8 mags, yields a distance 
~ 700 pc. If we assume a model with the same ([5.8 fim]- 
[8.0 /im]) color but which is slightly more inclined, with 
a [3.6 /im] ~ 9.7 mags, we obtain a distance ^ 330 pc. 
These distance estimates are highly speculative, as the 
uncertainty in the magnitudes makes the observed colors 
marginally consistent with later- type models. Further- 
more, these models are relatively untested. 

In this context, i.e., determining a distance to CB190, 
source 1 (see fig.S]) draws attention for two reasons. This 
star has a large amount of diffuse emission at 24 and 
70 /im. It is also coincident with the truncation of the 
^■^CO emission on the northern edge of CB190. Based on 
these two facts we conclude that source 1 is very likely to 
be associated with the cloud. Our spectrum of source 1 
shows it to be a B7 star, based on the strong Balmer ab- 
sorption lines and Hel features that are evident in fig. [H 
Assuming this is a main-sequence star, its distance is 
~ 400 pc, and if it is a giant (luminosity class III) its 
distance is ~ 600 pc. The broad feature observed in the 
spectrum of this source near 5700A may be due to the 
reddening curve. As a matter of historical interest, A. 
J. Cannon classified this star as B9. In the following 
analysis, where it is relevant, we assume a distance of 
400 pc. 

3.3. Extinction law analysis 

We have used the IRAC and 2MASS data to probe the 
extinction law in CB190. These data are used to generate 
the color-color plots shown in fig. [101 We calculate the 
errors in these colors by adding the respective IRAC and 
2MASS K-band errors in quadrature; the median val- 
ues of these errors are plotted in the lower right corner 
of fig. [ini To analyze the colors, we co mpare the best-fit 
redden ing vector to those measured bv[l ndebetouw et al.l 
()2005( ). We consider stars to be reddened if their colors 
are more than 2a away from the mean colors of the en- 
semble. These reddened stars are indicated in fig. [TUl as 
open boxes. We then fit the slope of the reddening vector, 
including both x- and y-er rors, using the IDL routine fi- 
texy.pro ( Press et al.lll99^ . We also include the median 
value of the unreddened stars and assign it zero error. 
The resulting reddening vector slope is shown in fig. [10] 
as a solid bla ck line. Because this slop e agrees reasonably 
well with the llndebetouw et al.l (|200El) extinction vectors, 
we conclude that the dust found in CB190 is not anoma- 
lous and that we are justified in using a "standard" red- 
dening law to determine the d ust properties in th i s glob - 
ule. The agreement with the llndebetouw et al.l (|2005l ) 
extinction is not s urpri sing because, as was sh own in 
'H arris et afl ([1978') and'Rieke fc Le bofskvl (|1985( ). there 
is generically almost no difference in the infrared extinc- 
tion in dense clouds even though the visual bands may 
deviate significantly from their low-density values. We 
note that we exclude source 2, discussed in § 3.1, from 
this analysis for two reasons: first, because it has a very 
large 24 //m excess, and second, because its broad-band 
SED is consistent with the colors of a proto-star, being 
too red to be simply due to a foreground dust screen and 
a normal star (see fig.lH]). 

3.4. The 24 fJ-m, shadow: Optical depth and column 
density profile 
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The 24 fim image shows a shadow, or depression in 
the emission, centered at RA = 19''20™48^ Dec = 
+23°29'45" (see fig. [J), coinciden t with the location of 
the dark cloud CB190 ( L771; e.g.. IClemens fc BarvainisI 
ll988HVisser et al.l l2002^. This shadow, about 70" in ra- 
dius, is coincident with the peak in the ^■^CO and ^^CO 
maps (see figs. [3] and S]), the 160 ^m emission (see fig.H]), 
and the SCUBA emission (see fig.|6]). We do not observe 
a shadow at 70 and 160 fim. In fact, at 70 /xm the cloud 
is at best only marginally detected in emission and may 
be dominated by the light from source 1 just to the north 
of the 24 /zm shadow. The shadow at 24 fj,m is a result 
of the cold and dense material in CB190 blocking the 
background radiation. 

In the following derivation of the density profile we cor- 
rect for the large-scale foreground emission at 24 (see 
foe and discussion below). This large-scale component 
is likely to be composed mainly of zodiacal light; if we 
underestimate it we will effectively wash out the shadow 
signal and therefore underestimate its optical depth pro- 
file, density profile, and mass. Conversely, if we over- 
estimate the large-scale background level, we will also 
overestimate the optical depth profile. Keeping this in 
mind, we estimate the zodiacal contribution by using the 
darkest parts of the image. Therefore, strictly speaking, 
we are deriving an upper limit to the optical depth pro- 
file. This approach is conservative because it allows us 
to derive a robust profile consistently and independently, 
without having to use other data to set the normaliza- 
tion of the density profile. While it is possible that the 
emission varies on shorter scales and, more specifically, 
emission from CB190 fills in the 24 fim shadow, we con- 
sider this to be unlikely due to the fact that the shadow 
is not detected at 70 /im. It is not plausible that CB190 
would be emitting significantly at 24 /im while remain- 
ing undetected at the longer wavelength. In the following 
discussion we describe the method used to derive the op- 
tical depth and extinction profile for the 24 fim shadow. 

First, we estimate the overall large-scale uniform back- 
ground level in the image, foe- We do so using two 
dark regions in the image, free from sources, indicated 
by boxes in fig. [Tl one to the North-East and one to 
the South- West of the shadow, each one about 50" on a 
side. We use the first percentile flux value in these boxes, 
fuc = —5.732 mJy arcsec"^, as a lower limit on the over- 
all uniform level in the image. We use this background 
level to set the true image zero level by subtracting it 
from the original image. Then we mask out all bright 
sources in the image by clipping all pixels with values 
3(7 above the mean. Finally, we use this background- 
subtracted and bright-source masked image to derive an 
optical-depth and extinction profile. For completeness, 
we estimate the error in fjjc by simulating the pixel dis- 
tribution in the two dark regions. We estimate anc by 
simulating Upi^ — 3362 pixel values drawn from a nor- 
mal density with mean and standard deviation equal to 
those measured in the two regions and storing the first 
percentile pixel value. We repeat 10^ times, and cal- 
culate the standard deviation in the simulated fjjc of 
CDC = 0.083 mJy arcsec"^. 

We proceed by dividing the shadow into regions of 
nested (concentric and adjacent) annuli 2.5" in width, 
the inner-most region being a circle with a radius of 2.5". 
These regions are centered on the darkest part of the 



shadow, which is also roughly coincident with the peak 
in the ""^^CO emission. We measure the average flux per 
pixel in each region out to a radius of ~ 100", chosen to 
be big enough to allow for a reasonable estimate of the 
background immediately adjacent to the shadow. We 
show the derived radial profile in fig. [11] Based on this 
profile, we set the boundary of the shadow to be at a 
radius of 67"5 (marked in fig. [11] with a solid line) where 
there is a flattening in the derived profile. We estimate 
the background flux by averaging the values of the annuli 
outside 70" and within 100". This background value is 
marked in figure [TT] with a dashed line. 

We use this profile to calculate the optical depth T24 
of the shadow. In a given annulus T24 is given by T24 = 
— ln(///o), where Iq is the background level and / is 
the shadow fiux. The calculated value of T24 varies from 
~ 1.5 at the center to about 0.02 at a radius of 67'.'5. 
The average mass column density in each annulus is then 
given by 

E = ^^/, (2) 

Kabs,2A 

where /(= 100) is the gas-to-dust ratio, and Ka(,s,24 is 
the absorption cross-section per mass of dust. In this 
work we use th e value of = 5.283 x 10^ cm^ gm^^ 

calculated by iDraind (2003a, b) from his Ry ~ 5.5 
model. The choice of a model with a high i?y value 
relative to the diffuse ISM is intended to account for 
some of the grain-growth effects likely to be taking place 
in CB190, as evidenced by the depletion indicated by 
our ^^CO and ^^CO data. This column density pro- 
file can be converted to an extinction profile using the 
relation N(H 2)/Ay = 1.87 x 10^^ atoms cm"^ mag~^ 
(Bohlin ct al. 1978), where we assume that all of the hy- 
drogen is in molecular form. We use a value of Ry = 3.1, 
appropriate for the diffuse ISM where this relation was 
measured, to convert from E{B — V) to Ay- The corre- 
sponding extinction is given by 



1.87x10^^ cm 2 mag ^/iH2'™H 



(3) 



where fin^ is the effective molecular weight per hydrogen 
molecule. The molecular weight per hydrogen molecule 
is related to the mass fraction of hydrogen, A^(II)/A^ 
(where 7W = 7W(H) -f 7W(He) + A1(Z)), by 



M 



2M 



2M 



mH^(H2) mYiN{B) X(H) 



(4) 



For a cosmic hydrogen mass fraction A^(II)/A^ = 0.71, 
/iHa = 2.8 (Jens Kauffmann - private communication, 
2006). We consider two sources of uncertainty in the 
column density profile for the 24 /im shadow: the DC- 
background error and the uncertainty in our assumed 
dust model, Kabs, 24. We measure an error in the DC- 
background level of aoc = 0.083 mJy arcsec^'^ (de- 
scribed above) . We then calculate the corresponding un- 
certainty in Kv by propagating udc through the equa- 
tion for T24 and scaling appropriately. We estimate the 
dust model uncert ainty to be half of the differ ence in 
Kafcs, 24 between the IWeingartner fc Draind (|2001[ ) Ry — 
3.1 model and the Ry = 5.5 model. To obtain an esti- 
mate of the total uncertainty in Ay, we add these two 
components and a 10% systematic error floor in quadra- 
ture. 
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The column density profile and corresponding errors 
are shown in figure [12] along with individual 2MASS 
point-source extinction estimates. We use 2MASS 
sources with data quality flags better than, and includ- 
ing, "UBB" , and include "C" quality measurements in 
a bandpass when the other two fllters are "B" quality 
or better. Sources with upper limits in two or more 
bands are rejected. This selection ensures that we in- 
clude the maximum number of quality extinction mea- 
surements while not biasing the object selection towards 
stars with lower amounts of reddening. The reddening 
of the accepte d sources is measu r ed us ing the [ J — K) 
color and the iRieke fc Lebofskvl (|1985| ) extinction law: 
Av = E{J -K)/Q.ll. The "intrinsic" (J - K)o = 1.34 
color of the stars is measured in the 70" to 100 " an- 
nulus centered on the shadow, to isolate the effects of 
the globule from that of nearby diffuse dust. We reject 
sources between 0" and 70" with E{J — K) values lower 
than —0.37, the l-tr value for the scatter in the sources 
between 70" and 100". This l-cr range is marked in fig- 
ure [13 with two dotted lines. Two sources fulfill this 
criterion, located at r ~25", 67" from the center. We re- 
ject one other source with an Ay ~ 3 located at r ~ 9" 
on the basis that it is likely to be a foreground object. 
We calculate a probability of 1.4% of finding three fore- 
ground stars, using the stellar nu mber densities from the 
Nearby Stars Database (NStars: iHenrv et al.l[2003^ . It 
is not surprising that this probability is so low because 
the NStars catalog is not complete. A more accurate es- 
timate of the likelihood of finding three foreground stars 
within the shadow is determined as follows. In the 70" 
to 100 " annulus there are three stars with upper limits 
below Av = and two stars with Ay < —0.37. Further- 
more, the 70" shadow and the surrounding 70" to 100 " 
annulus have the same area. We therefore expect ~ 4 
foreground stars in the shadow region, and hence reject- 
ing three stars is reasonable. Using the remaining stars 
we calculate the best-fit Gaussian parameters for the dis- 
tribution of Av values in two annuli, from " to 40 " 
and from 40 " to 70". Lower limits are not treated dif- 
ferently than proper detections in our fitting procedure. 
Therefore, when we plot the mean values from these fits 
versus average radius the values are shown as lower lim- 
its (see fig. [T2|). The trend of increasing column, or Ay, 
with decreasing radius can be seen clearly. 

In this section we have presented a new technique for 
analyzing 24 fim shadows which allows for a smooth esti- 
mate of the density profile of these cold cloud cores at a 
6" resolution. Most importantly, this method traces gas 
and dust down to the densest regions in the cloud. 

4. MASS ESTIMATES 

In the following mass estimates we assume a distance 
to CB190 of 400 pc (cf. § 3.2). We present a summary 
of these calculations in Table 1. 

4.1. 24 IJ-rn shadow mass 

We use the optical depth profile derived in § 3.4 to 
calculate the mass of the 24 /im shadow. The dust mass 
in a given annulus is 

Md = — ^^rip„np„i:»^, (5) 

Kabs,24 

where T24 is the optical depth in the annulus, Kafcs, 24 
is the absorption cross section per mass of dust at 



24/im, flpix is the solid angle subtended by a pixel, 
Hpix is the number of pixels in the annulus and D is 
the distance to the cl o ud. U sing the Kabs,24 from the 
IWeingartner fc Draind ()2001h Milky Way synthetic ex- 
tinction curve model with Ry = 5.5, a gas-to-dust ratio 
/ = 100, and summing over all the annuli within 70", 
we derive a total 24 /xm shadow mass of 9.7Mo (see in- 
set of fig. [T2)) . Using equation 4 and standard propa- 
gation of errors we find that the fractional error in the 
mass is equal to 3.3 times the fractional error in the lo- 
cal background flux level, Iq (marked as a dashed line in 
flg. [TT|) . For a fractional error in the background level 
of 4%, corresponding to the the 95% confidence interval 
for the distribution plotted in fig. [TTJ we calculate a 14% 
error in the mass. We note that this error analysis does 
not include other significant systematic sources of error, 
like differences in calculated dust opacities. For example, 
lOssenkopf fc Henningj (|1994[ ) calculate a 24 /xm Kabs, 24 = 
8.69 X 10^ cm^ gm~^ (for the thin ice mantle model gen- 
erated at a density of n = 10^ cm~^) , a fac tor ^ 1.6 
greater than the .Weingartner fc Draind (|2001h opacities 
used here, which would reduce our mass to ~6 M©. 

We adopt a mass of ~ 10 Mq, although the estimates 
from the 160 /im data imply this value may be a lower 
limit. This mass is ~ 10 times bigger than the mass de- 
rived from the CO observations (see § 5.1). Furthermore, 
for a spherical cloud of the same size and mass as CB190, 
and a temperature of 10 K (consistent with both the CO 
line data and the far infrared to submillimeter contin- 
uum fit) we derive a Jeans mass ~4.1 Mq. While such 
a low mass would imply instability to collapse if we only 
consider thermal pressure, this analysis neglects alter- 
nate forms of support - namely turbulence and magnetic 
fields - which likely play a significant role in CB190. 

4.2. 160 ijm mass estimate 

Assuming that the dust is optically thin, the dust mass 
associated with the 160 /im emission is given by 



B,^{T)Kabs,160 ' 

where D is the distance to the cloud and Kofcs.ieo is the 
absorption cross section per mass of dust at 160/im. We 
use a value for Kabs.iea = 40.14 cm^ gm~^, given by the 
Ossenkopf fc Henning ( 1994) model with thin ice mantles 
and generated at a density of n = 10^ cm~'^. These mod- 
els account for grain growth effects, which are likely to 
be taking place inside the cold, dense cloud environment. 
These authors show that the effects of grain growth on 
the opacities are, however, minor. To calculate the dust 
mass we integrate the 160 /im flux within a 6.25 pixel 
(100") radius aperture centered on the 24 /j,m shadow 
(and the peak ^^CO emission) and subtract an estimated 
pixel background level. We estimate this background per 
pixel by averaging the flux in an annulus with an inner 
radius of 20 pixels and a width of 4 pixels. We assume a 
gas-to-dust ratio / = 100. A typical temperature for dust 
grains heated b y the i nterstellar radiation field (IRF) is 
-18 K Id raindl2003aD : if the cloud is externally heated 
by the IRF then the total mass, given by equation 4, 
is ~ 1.2 Mq. However, we have shown in § 2.5 that 
the temperature of the cloud is in the range of 10. 4K to 
12K, giving a total mass estimate between ^--^45 Mq and 



7 



--14 Mq. The lOssenkopf fc Hennind ()1994f ) models also 
include opacities for grains with thick ice mantles; this 
opacity would sc ale these masses down by a fa ctor of 1.3. 
We note that the lWeingartner fc Drai^^ (|200H ) Ry = 5.5 
ISM opacity of Kabs,i60 = 10.12 cm^ gm~^ yields masses 
that are larger by a factor of 4. However, these models 
do not take into account the grain growth effects that are 
thought to occur in cold, dense regions. We note that the 
temperature of CB190 is not likely be constant through- 
out the cloud but instead probably decreases inwards; 
this gradient would bias the 160/im observations toward 
hotter dust on the outside of the globule and hence we 
might underestimate the mass. 

4.3. CO mass estimate 

We use observations of ^^CO and ^'^CO J = 2 ^ 1 to 
estimate a mass for CB190. We use the ^^CO data to es- 
timate an excitation temperature in the cloud; this tem- 
perature is then used in conjunction with the ^^CO data 
to estimate a mass, using standard assumptions (e.g., 
[Rohlfs fc W ilson 2004). We derive a ^^CO gas tempera- 
ture ~ 10 K. For an isotope ratio of [i^CO/^CO] = 50 
and a stan dard non-depleted r atio of n{CO)/n{H2) — 
1.5 X 10-4 (iKulesa et al.l[2?)05h we derive a cloud mass 
M = 4.2 Mq by summing over the entire ~ 10' by 10' 
area of the image. For an aperture with a radius of 70" 
centered on the peak of the ^^CO emission (and the 24/im 
shadow) we derive a total gas mass of ~ 1.1 Mq. This 
CO mass is much less than the 9.7 Mq derived from 
the 24/im shadow profile (cf. § 4.1). These contradic- 
tory mass measurements can be reconciled if freezeout 
is an important effect, which we consider to be a likely 
scenario, or if the n{C0)/n{H2) ratio is over-estimated 
by about an order of magnitude. We note that masses 
derived in this fashion will simply scale with the isotope 
ratio. For comparison we calculate the virial mass of the 
globule. Using our higher resolution data we measure the 
FWHM of the i^CO line to be 0.97 Km s'^ at the lo- 
cation of the shadow. However, the C^*0 data suggests 
that there may be two velocity components broadening 
the "CO hue. The hue width of the C^^O observation 
is roughly ISV ^ 0.5 km s"^ at FWHM, in agreement 
with the NH3 linewidth. However, like the ^^CO line, 
it appears to have some asymmetry on the red side (see 
fig. E]). If we estimate the total width of the line by dou- 
bling the HWHM from the blue side, then we obtain a 
/SV ~ 0.2 km s^^. We take 0.5 km as an upper 
limit on the line width and 0.2 km as a more plausi- 
ble value. These line widths correspond to a virial mass 
between IS.OMq and 2.9Mq, respectively. 

5. BONNER-EBERT MODELS AND POSSIBLE SUPPORT 
MECHANISMS 

When analyzing cold cloud density or extinction pro- 
files it is common to use theoretical density profiles 
to provide physical insight to the system(s). Bonnor- 
Ebert models are one such choice; they are solutions 
to the isothermal equa tion, also known as the m odified 
Lane- Emden equation ()Bonnoilll95l lEbertlll955h . The 
isothermal equation. 



coordinate, where 



(7) 



describes a self-gravitating isothermal sphere in hydro- 
static equilibrium. Here, — r/vc'is the scale- free radial 



kT 



{jAirGfirnHPc 



1/2 



(8) 



is the scale-radius, r is the physical radial coordinate, 
G is the gravitational constant, fj, (= 2.37) is the mean 
molecular weight per free particle, ruH is the mass of a 
hydrogen atom, k is the Boltzmann constant, T is the 
temperature, and pc is the central mass density. Finally, 



* = -ln(p/pc), 



(9) 



where p is mass density as a function of radius. Equation 
7 can be solved, with appropriate boundary conditions, 
to obtain the scale-free radial profile of a Bonnor-Ebert 
cloud. The singular isothermal s phere (n (x r~^) is a 
limiting solution, with S^max 00 (|Chandrasekhajl939r ). 

Any given solution to equation (7) is characterized by 
three parameters: the temperature T, the central den- 
sity Pc, and the outer radius of the cloud R. Once an 
outer radius is specified, a model will be in an unstable 
equilibrium if ^max > 6.5, where 
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R/rc 



(10) 



Equivalently, an unstable model will have a density con- 
trast betwe en the center and the edge of the cloud greater 
than 14.3 (jBallesteros-Paredes et all I2003D . For_ more 
discussion of Bonnor-Ebert models see, e.g.. lEvans et al 



(20011 iBallesteros-Paredes et al.1 (I2003D. iHarvey et al 
(j2003E . lLada et al.l (|2004l ^ and lShirlev et all (|2005f ). 



5.1. Bonnor-Ebert fits 

We generate a set of Bonnor-Ebert models over a 3-D 
grid in temperature, outer radius, and central density. 
Our temperature grid ranges from T = 8.0 K to 98.0 K 
in steps AT = 2.0 K. Our radial grid varies from Omax = 
45'.'0 to 95'.'0 in steps of AO^ax = 2'.'0, and we assume a 
distance of 400 pc. For the central density grid, we vary 
the ^max parameter from £,max — 3.0 to 19.0, in steps of 
AS^max — 0.2. As can be seen from equations (8) and 
(10), at fixed temperature and 9max {— R/D, where D 
= 400 pc is fixed), varying S^max is equivalent to varying 
the central density pc- We then integrate our calculated 
density profiles to obtain column density profiles: 



NBE{r) = 2x / p(r') 



dr' 



(11) 



where r is the projected distance from the center of the 
shadow and the integration is along the line of sight r' . 

We calculate the best-fit model by finding the mini- 
mum over the grid in temperature, Omax, and £,max, 
where 

{N, - NsE^i)'' 



X 



(12) 



Here, iV^ is the measured column density, ai is the corre- 
sponding uncertainty in Ni, NBE,i is the Bonnor-Ebert 
model extinction, and the sum over i represents the inte- 
gration over the spatial coordinate. Errors in the best-fit 
parameters are calculated using a Monte Carlo approach; 
we generate a set of 10"^ Ay extinction profiles, given by 
equation 3, with values of Jbc and Kabs, 24 with normal 
distributions, described in § 3.4. We also include the 10% 
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systematic noise floor in our calculation of the mock Ay 
profiles. For each one of these artificial column density 
profiles we fit a model according to the procedure out- 
lined above. The errors in the fitted parameters that we 
quote are the standard deviations in the resulting best-fit 
parameter distributions. 

We find a best-fit Bonnor-Ebert model with a tem- 
perature T = 56.0 ± 7.3 K, Bmax = 69'.'0 ± 0.1, and 
£,max — 7.2 ± 0.2. The central density of this model 
is Tic = 1-78 X 10^ cm~'^. We show this model with the 
data in fig. [131 Despite the reasonably good agreement 
between the model and the 24 /xm profile, this fit is in- 
consistent with the molecular data and the SCUBA data: 
the temperature from the CO line is ~ 10 K while the 
continmmr fit yields a temperature of 10 to 12 K. We 
note that reducing the assumed distance to CB190 low- 
ers the best-fit model temperature: a distance of 200 pc 
brings the model temperature down to ~28 K. However, 
we rule out this short distance based on the likely asso- 
ciation of star 1 with CB190, see fig. ID We discuss this 
distance in more detail in § 3.2. Even though the best-fit 
£,max is only shghtly greater than 6.5, our best-fit mod- 
els always have either inconsistently high temperatures 
or unacceptably small distances. Therefore we must rule 
out stable Bonnor-Ebert profiles as accurate representa- 
tions of CB190. 

5.2. Turbulence 

The temperature of CB190 is ~10 K, based on our 
160 A^ m, 850 jum, and CO observations, typ ical for cold 
cores (jLemme et al.lll996HHotzel et al.ll2002[ ). This tem- 
perature corresponds to a thermal line width of AVth ^ 
0.09 km s~^ for C'^^O. As discussed in § 4.3, the observed 
line width of the C^^O observation is between 0.2 km s^^ 
and 0.5 km s~^, much broader than that expected from 
thermal support. We take 0.5 km s~^ as an upper limit 
on the turbulent line width and 0.2 km s^^ as a more 
plausible value, 

From Hotzcl ^t al.l (|2002[ ). the energy contributed by 
turbulence to the support of the cloud is 



30 VAT/, 



2 - Eth 



th 



(13) 



Assuming that /W = {0.5,0.2} km s"^ then 
Eturb/Eth — {3.5,0.5}. In this case, the maximum 
amount of energy is Etot = ^-^Et, which would supply 
support equivalent to a temperature of 45K while the 
more plausible value is only ^ 10 — 15 K. These temper- 
atures are lower than our Bonnor-Ebert best-fit value of 
T — 56.0 ± 7.3 K. Although turbulence might in fact be 
a significant source of support in CB190, it is not likely 
to provide enough outward pressure to prevent collapse. 

5.3. Magnetic support 

Here we consider the effects of magnetic fields, an alter- 
native to turbulence as a support mechanism in CB190. 
Under a b road range of condi tions, the magne t ic pre ssure 
is B^/%TT (lBosslll997D . From lStahler h Fallal ((200l . the 
mass that can be supported given a cloud radius R and 
magnetic field magnitude B is 



CB190 has an average column density of {N) = 2 x 
10^^ cm^^, which corresponds to a line-of-sight magnetic 
field of Bios = 50 ^G, given by the observed correlation 
for 17 cloud s with confirmed magnetic field detections 
(jBasul l2004[ ) . We note that this relation does have a 
large scatter of about 0.2 dex. Assuming equipartition, 
the total magnitude of the expected magnetic field is 
then Btot = x Bios — 87 fiG. Using a radius of 
R ~ 0.13 pc, equation 13 gives a mass M = 11 Mq. 
Therefore, the magnetic field could be of sufficient mag- 
nitude that it may support this cloud and retard collapse. 
This result is interesting because magnetic support is of- 
ten overlooked in studies of cold cloud cores, yet in the 
case of CB190 it may play a dominant role. We note that 
a decrease of 0.2 dcx in the magnetic field will make it 
insufficient to support CB190. 

5.4. Comparison with other globules 

iKandori et al.l ()2005l ) summarize observations of the 
density structure of dark globules. They show that, 
of 11 starless cores with good density measurements, 
7 have profiles consis tent with purely thermal support. 
iTeixeira et al.l (|2005f ) report three additional starless 
cores in Lupus 3, of which two appear to be supported 
thermally. That is, of 14 such cores, 9 have profiles con- 
sistent with pure thermal support. In general, the five 
cores where an additional support mechanism is required 
do not have adequate line width measurements to assess 
the role of turbulence. 

Our study of CB190 places it among the relatively rare 
class of cores that cannot be supported purely thermally. 
Our high resolution line measurements indicate that the 
turbulence in this core is inadequate for support also. It 
is plausible that the magnetic field in the globule sup- 
plies the deficiency, if it is at an average level measured 
for other cold cloud cores. If CB190 is currently sup- 
ported in this way, it is at an interesting phase in its 
evolution. For the properties of this clo ud, the ambipo- 
lar diffusion timescale is ^ 3 x 10^ yr (IStahler fc Fallal 
1200 5'). which is about a factor of 10 longer than the free- 
fall collapse timescale. Thus, it is predicted that magnet- 
ically supported cores are unstable over about ten mil- 
lion years, as their ma gnetic fields leak out t hrough am- 
bipolar diffusion (e.g., iCrutcher et "all 119941; B oss 1993; 
Indebetouw & Zweibel '20001; ISigalotti fc KlappI I2000l ; 
Tassis fc Mouschoyi as 2004). At that point, we would 
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expect CB190 to begin collapsing into a star. If the mag- 
netic field is lower than average, this process may already 
have begun. 

6. SUMMARY 

We have combined Spitzer MIPS and IRAC data with 
HHT and GBT millimeter data of CB190 and arrive at 
the following conclusions: 

• We introduce a method for studying the structure of 
cold cloud cores from the extinction shadows they cast 
at 24 /im. 

• We derive an Ay profile of the 24 ^m shadow that is 
in good agreement with the reddening estimates derived 
from the 2MASS data at the outer edges and reaches a 
maximum value of ~ 32 visual magnitudes through the 
center. 

• The mass measured from the optical depth profile is a 
factor of ^--^2 greater than the Jeans mass for this object. 
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• We fit Bonnor-Ebert spheres to our Ay profile and find 
tliat the best-fit temperatures are in contradiction with 
the CO observations and the thermal continuum data, 
which indicate much lower temperatures for this globule. 
We also show that turbulence is probably inadequate to 
support the cloud. However, magnetic support may be 
enough to prevent collapse. 

These pieces of evidence together form a consistent pic- 
ture in which CB190 is a cold dark starless core. Al- 
though collapse cannot be halted with thermal and tur- 
bulent support alone, the magnetic field may contribute 
enough energy that it could support CB190 against col- 
lapse. Hence, magnetic field support should be in- 
cluded in evaluating the stability of other cold cloud 
cores. CB190 appears to be at an interesting evolution- 
ary phase. It may be in the first stages of collapse (if 
the magnetic field is weaker than average). Alternately, 
if it is currently supported by magnetic pressure, it is 
expected that collapse may begin in some ten million 
years as the magnetic field leaks out of the globule by 
ambipolar diffusion. 
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TABLE 1 

Summary of mass estimates 



Data 


Temp. [K] 


Aperture size ["] 


Mass [Mq] 


24/im 




70 


10 


Jeans Mass 


10 


70 


4 


160/im" 


10 


100 


~ 45 




12 


100 


~ 14 


"CO & i^CO 




70 


^ 1* 


C^*0 virial mass'' 




70 


3 



Evidence for CO freezout implies that this value is a lower limit. 
^ Indicated as upper limits to allow for grain growth in dense regions. 
^ Based on an estimated line width of 0.2 km s~^. 




Fig. 1.— Gallery of CB190 data. The images arc centered on the 24/im shadow at RA = 19''20'"48^ Dec = +23°29'45", and are 
oriented such that north is up and east is to the left. The images have a height of 24' and a width of 12'. The corresponding wavelengths 
are labeled in the top right corners. The optical image is from the red Digital Sky Survey. The original MIPS mosaic pixel scales arc 
~ 1^24 pix^-', 9''8 pix^-'^, and 16" pix^^, for 24 /im, 70 /im, and 160 /im respectively. The 24 /im image displayed here is binned down by 
a factor of 4, and the DSS, 70, and 160 /im images arc re-gridded to the same scale. All four images are marked with a 100" circle centered 
on the location of CB190. The 24 fim image is marked with two ~ 50" X 50" boxes indicating the background regions used to estimate the 
overall uniform background level (see § 3.3 for more details). The 160 //m image shows the aperture used for the dust mass determination 
(see § 3.2) indicated a black circle with a radius of 100", and the inner and outer sky annuli radii at 5.'3 and 6.'4, both marked as white 
dashed circles. 
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Fig. 2.— Gallery of IRAC CB190 images centered on the 24/im shadow: RA = 19''20'"48'', Dec = +23°29'45". The images are oriented 
such that north is up and east is to the left, they are 6' on a side, and the corresponding effective wavelengths of the four IRAC channels 
are labeled. The original mosaiced IRAC image pixel scale is 0'.'6 pix~^ and the images displayed here are shown at 1'.'2 pix~^. The black 
circles are centered on the 24^m shadow and have a radius of 100". 




-4-2024 -4-2024 
Aa (arcmin) Aa (orcmin) 



Fig. 3.— Maps of the integrated ^^CQ (J = 2-1, v = 220.399 GHz) and ^^CO (J = 2-1, v = 230.538 GHz) lines taken with a FWHM 
telescope resolution of 32" and convolved to a square grid cell spacing of 16". The central coordinates of the images are RA = 19'"20™49^, 
Dec = -|-23°29'57". The color-bars indicate the temperature scales in each map. The ^^CO contour levels are 2, 3.5, 5, and 6 K, and the 
i^CO contour levels are 0.8, 1.4, and 2.0 K. 
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Aa (arcmin) 



Fig. 4. — 24 fim image with "'^^CO contours overlayed. The image is ~ 10' on a side, is centered on the CO map at RA = 19'"20'"48=, 
Dec = +23°29'57", and is displayed at the original mosaic pixel scale of 1'.'24. The ^'^CO contour levels are 0.8, 1.4, and 2.0 K. The white 
circle is centered on the 24^m shadow and the peak of the ^^CO emission. The two arrows indicate two sources that are likely associated 
with CB190; their broad-band SEDs are shown in fig. |8] and discussed in § 3.1. 
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Fig. 5. — High resolution spectra of ^^CO and ^'^CO J = 2 - 1 lines at the position of the peak CO optical depth, which coincides with 
the center of the 24 lira shadow. The spectral resolution of the ^"^CO line is 0.25 MHz or 0.34 km s~^, and the resolution of the ^^CO 
line is 0.05 MHz or 0.061 km s~^ (R = 5,000,000). The vertical scale of the ^^CO spectrum is multiplied by 2.5. Total integration time 
was 15 min on source. RMS noise in the ^^CO spectrum is 0.12 K (main-beam brightness temperature. The dotted line shows the C'^^O 
J = 2 - 1 spectrum obtained with the Caltech Submillimeter Observatory at a spectral resolution of 0.1 MHz or 0.13 km s"'^ (Michael M. 
Dunham - private communication, 2006). The vertical scale of the C^*0 spectrum is multiplied by 6. 
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Fig. 6. — 160/im image detail, centered on the 24^m shadow, with 850/im contours overlaycd in black. Tiic 160/im image is shown at the 
original mosaic pixel scale of 16'.'0 pix"'^. The 850/im contour levels are 0.1, 0.2, 0.3 and 0.4 Jy beam-\ and the beam is a 32" FWHM 
Gaussian. The grey circle indicates the location of the 24 /im shadow and is 48" in radius. 




Fig. 7. — 70/im, 160/im, and 850/im photometry using a 48" radius aperture centered on the 24^m shadow coordinates. We show two 
models for the dust emission with the indicated assumed values of (3 and the corresponding best-fit temperatures. 
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Fig. 8.— vs. A for three stars: star 1 (HD344204, top), HD1608 (middle), and Vega (bottom). The data extend from 3615A to 6900A 
at 9A resolution. 
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Fig. 9. — Broad-band SEDs of two sources associated with CB190; wc show the 2MASS, IRAC and MIPS magnitudes. The x's show the 
magnitudes of star 1 (HD344204), at RA = 19''20'"47^ Dec = 23°31'40.637". The solid circle indicates the magnitude of star 1 using a 
large aperture to include the dust emission surrounding the source. The open squares show the magnitudes of star 2, at RA = 19'' 20^157", 
Dec = 23°31'37.6". 
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Fig. 10. — Color-color plot of stars observed towards CB190. We combine 2MASS K-band data with 3.6 fim and 5.8 IRAC channels. 
The squares indicate the reddened stars, selected with a 2-a clipping criterion, and the filled circle s indica t e the unreddened stars. The 
black arrow shows the best-fit extinction vector slope of 1.2, in good agreement with the Indcbctou w et aLl (120051) extinction slopes. The 
median errors for the data are shown in the lo wer right. The extinction vector derived from the K-[3.6/im] vs. K-[4.5/tm] colors also shows 
good agreement with llndebetouw et al.l l|2005l) . 
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Fig. 11. — The radial flux profile for the 24 fim shadow. The squares indicate the mean counts per pixel in each 2''5 annulus. The chosen 
truncation radius for the shadow, at 67''5 is marked as a solid line; the local background flux level /q is marked with a dashed line. 
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Fig. 12. — Compeirison of the CB190 column density profiles: the 24 /Lim shadow profile is shown with dark grey squares, and the 
individual 2MASS stellar estimates are marked with light grey filled circles. The two dotted fines indicate the l-cr scatter in the 
2MASS data just outside the globule, between 70" and 100". The black symbols indicate the best-fit Gaussian mean values for the 2MASS 
data in three radial bins, inside 40", 40" to 70", and 70" to 100", and are plotted versus the average radius in each respective bin; the 
two inner-most points are shown as lower limits (bla<;k arrows) because of the inclusion of individual 2MASS lower limits in the best-fit 
Gaussian calculation. The inset shows the corresponding enclosed mass for the 24 /nm extinction profile. 
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Fig. 13. — The 24 urn column density profile (squares) is shown with the best-fit Bonnor-Ebert model (solid line). The errors in the data 
are about the size of the squares. The best-fit model parameters are indicated: temperature T = 56.0 K, dmax = 69"0, and ^max = 7.2. 
The right-hand-side vertical axis indicates the magnitudes of visual extinction, or equivalently, T24 , the optical depth at 24 fim. 



